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When the constituent spins have an energetic preference to lie along an easy-axis, triangular 
and Kagome lattice antiferromagnets often develop long-range order that distinguishes the three 
sublattices of the underlying triangular Bravais lattice. In zero magnetic field, this three-sublattice 
order melts either in a two-step manner, i.e. via an intermediate phase with power-law three- 
sublattice order controlled by a temperature dependent exponent r](T) G (|, |), or via a transition 
in the three-state Potts universality class. Here, I predict that the uniform susceptibility to a small 
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easy-axis held B diverges as y(i?) ~ \B\ 4-9r ' in a large part of the intermediate power-law ordered 

phase (corresponding to r/(T) G (|, |)), providing an easy-to-measure thermodynamic signature of 
two-step melting. I also show that these two melting scenarios can be generically connected via an 
intervening multicritical point, and obtain numerical estimates of multicritical exponents. 

PACS numbers: 75.10.Jm 


In frustrated antiferromagnets mm, magnetic ions 
(spins) form a lattice whose geometry causes the dom¬ 
inant antiferromagnetic interactions between neighbours 
to compete with each other. This allows weaker further- 
neighbour interactions or quantum fluctuations to se¬ 
lect complex patterns of spin order at low temperature. 
Models of frustrated easy-axis antiferromagnets H®, in 
which spins can lower energy by orienting along a fixed 
axis, provide interesting examples of this behaviour. 

Such models are also relevant in other experimental 
contexts. For instance, the low-temperature behaviour 
of monolayers of adsorbed gases on substrates with tri¬ 
angular symmetry mm lias been modeled |T2] in terms 
of a triangular lattice of Ising spins a z s = ±1 (z com- 

ponents of spin-half moments S^ = erg/2) with antifer¬ 
romagnetic Ising interactions Ju%a z s between nearest- 
neighbours [13] [14] and weak ferromagnetic Ising inter¬ 
actions between further-neighbours. More recently, the 
magnetic properties of honeycomb networks [MIS] of 
magnetic wires (dubbed artificial Kagome-ice) have been 
analyzed UMB in terms of a similar Ising model on the 
Kagome lattice [22] • In both examples, further-neighbour 
couplings cause the Ising spins to develop ferrimagnetic 
three-sublattice order at low temperature, i.e., freeze into 
a pattern which distinguishes the three sublattices of the 
underlying triangular Bravais lattice and gives rise to a 
small net moment along the easy axis. 

Several other easy-axis spin systems on triangular and 
Kagome lattices exhibit ferrimagnetic three-sublattice or¬ 
der [23lf34] . or closely related antiferromagnetic (no net 
easy-axis moment) three-sublattice order [35] , In zero 
field (B = 0) along the easy-axis, a Ginzburg-Landau 
theory E5H55] for the three-sublattice order parameter 
predicts that this ordering transition is described by a 
six-fold anisotropic effective model of ferromagnetically 
coupled XY spins [39], or, equivalently, by a generalized 
six-state clock model [40IHB1 . Rather unusually, such six- 


state clock models have multiple generic possibilities for 
continuous transitions: Order is lost either v ia a two-step 
melting transition, with an intermediate phase character¬ 
ized by power-law order [32], or via a sequence of two dis¬ 
tinct transitions, one of which is in the three-state Potts 
universality class and the other in the Ising universality 
class m mi- Perhaps motivated by this, the melting 
of three-sublattice order has been studied in a variety 
of triangular and Kagome lattice systems for over three 
decades now. In some examples [H [23H2H Eg, three- 
sublattice order is known to melt in a two-step manner, 
via a sizeable intermediate phase with power-law three- 
sublattice order controlled by a temperature-dependent 
exponent rj(T ) G (|, |). In other examples with fer¬ 
rimagnetic three-sublattice order, this order is lost via 
a three-state Potts transition, while residual ferromag¬ 
netism is lost via an Ising transition |19H21| . 

In this Letter, I analyze the melting of three-sublattice 
order in easy-axis antiferromagnets on triangular and 
Kagome lattices using a new coarse-grained description 
that explicitly keeps track of the uniform easy-axis mag¬ 
netization mode whose fluctuations are coupled to fluctu¬ 
ations of the three-sublattice order parameter. Using this 
description, which goes beyond the standard Ginzburg- 
Landau theory, I demonstrate that these two very dif¬ 
ferent melting processes can be generically connected via 
an intervening multicritical point M (Fig. [2]) with cen¬ 
tral charge [H] c^t G (1, §). Although the generalized 
six-state clock model correctly captures other generic 
ways [40] ]44] in which these two very different melting 
processes can be separated from each other in the phase 
diagram of such three-sublattice ordered systems, it fails 
to account for the existence of 3VL This underscores the 
importance of treating the uniform magnetization mode 
on the same footing as the three-sublattice order param¬ 
eter. 

I obtain numerical estimates of multicritical exponents, 
and argue that such multicritical melting may be experi- 
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mentally accessible in artificial Kagome-ice systems if the 
strengths of nearest and next-nearest exchange interac¬ 
tions can be increased relative to the long-range dipo¬ 
lar interactions. Additionally, for r]{T) £ (|,|) in the 
power-law ordered phase associated with two-step melt¬ 
ing, I show that the uniform susceptibility to a small easy- 
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axis field B diverges as x(B) ~ \B\ . I also argue 

that this easy-to-measure thermodynamic signature is of 
potential experimental relevance in the context of three- 
sublattice ordering of nearly-half-filled monolayers of ad¬ 
sorbed gases on triangular substrates, and in the context 
of experimental realizations of three-sublattice order in 
5 = 1 Heisenberg antiferromagnets with strong single-ion 
anisotropy on the triangular lattice. 

Order parameters and coarse-graining: I use the 
convention of Fig. 0 for labeling the sites [unit-cells] 
R = me x + ne y of the triangular [Kagome] lattice, and 
for labeling the three basis sites a = 0,1, 2 in each unit¬ 
cell of the Kagome lattice. With this convention, the 
complex three-sublattice order parameter if = |^>|e* e and 
the ferromagnetic order parameter M z are defined as: 
if = — e* = 5 L ( m + ri )5^ and M z = on the tri¬ 
angular lattice, while if = — „ e l? s L ^ m+n ~ a ' > S% and 

ii,Ot 

M z = J2 r a ^ Z R a 0n the Kagome lattice. Our coarse¬ 
grained description will be written in terms of an effec¬ 
tive Hamiltonian defined on a lattice whose sites r repre¬ 
sent clusters of spins of the original triangular or Kagome 
magnet. In this description, each cluster is characterized 
by an Ising variable Tp = ±1 representing the direction 
of the local easy-axis magnetization M c z luster , and by an 
angle Op that represents the phase of the local three- 
sublattice order parameter ^cluster- Comparison with 
long-wavelength properties of specific microscopic models 
is facilitated by choosing clusters that themselves form a 
coarse-grained triangular lattice, since this preserves the 
symmetries of the underlying triangular Bravais lattice 
in both triangular and Kagome lattice systems. 

Ginzburg-Landau theory: Let us begin by summa¬ 
rizing in this language the standard Ginzburg-Landau 
theory for three-sublattice ordering [5HIH58] : Transfor¬ 
mation properties of if under global spin-flip and lattice 
symmetry operations fix the form of the effective Hamil¬ 
tonian H xy for 9p. Leaving out certain chiral perturba¬ 
tions [MH3 that are not expected to be relevant 0S] for 
the transitions of the lattice magnets studied here, H xy 
may be written as 

H xy = — J xy cos (6p — 9p ) — h 6 ^ cos(66V) . (1) 

(rr') r 

where (rr f ) are nearest-neighbour links of our coarse¬ 
grained triangular lattice. The effective stiffness J xy > 0 
(encoding the energetic preference for three-sublattice 
order) and the six-fold anisotropy he, whose sign se¬ 
lects between ferrimagnetic three-sublattice order (with 


0 m = 2irm/6) and antiferromagnetic three-sublattice or¬ 
der (0 rn = (2 m + 1)7t/ 6), are both set by quantum fluc¬ 
tuations and subdominant further-neighbour couplings 
in the microscopic Hamiltonian. In this approach, the 
relative values of J xy and its higher harmonics (co¬ 
efficients of — cos {pOp — pOp ) for p = 2, 3) determine the 
nature of the melting process. These higher harmonics 
are omitted from H xy displayed above since they are not 
crucial for our subsequent discussion. 

New effective Hamiltonian: Next, I note that this 
standard Ginzburg-Landau description does not take into 
account the uniform magnetization mode whose fluctua¬ 
tions are coupled in a crucial way to fluctuations of the 
three-sublattice order parameter. This key observation 
leads me to a new coarse-grained effective model: 

H e ff // X y T -Rising Jqt ~ 7”y COS (3fly) , 

r 

_ (2) 

where Lfi s i ng = -Rising Tprp - h ^ t? , 

(rr') r 

with h oc B. To understand the rationale for the form of 
this effective Hamiltonian, it is useful to first note that 
L/eff has the same S3 x Z 2 symmetry as H x y , and re¬ 
duces, in the double limit he, Jg T —» 00 , to a generalized 
six-state clock model studied earlier 001 HU. However, 
the space of states at each site of H e g is enlarged by 
the presence of Tp to correctly account for the fact that 
the direction of Af z luster is correlated with the phase of 
^cluster, but not completely tied to it. The microscopic 
origin of various terms can now be understood as follows: 
Rising > 0 encodes the effect of subleading ferromagnetic 
interactions of the microscopic magnet, which tend to 
favour ferrimagnetic three-sublattice order. If he > 0, it 
is likely to be accompanied by a sizeable positive value of 
Rising in Hoff (since ferrimagnetic three-sublattice order 
corresponds to he > 0 in H xy ). Conversely, negative he, 
favoured by quantum-fluctuations in some examples 05], 
is likely to be accompanied by negligibly small Ji s i ng . The 
coupling Jg T > 0 correctly captures the fact that the 
values 9 = 0, 27r/3, 47t/3 (7t/3, tt, 5n/3), characteristic of 
ferrimagnetic three-sublattice order, are associated with 
a positive (negative) easy-axis magnetization, while the 
phase choices 9 = (2m + 1)7t/6, characteristic of anti¬ 
ferromagnetic three-sublattice order, are not associated 
with any net easy-axis magnetization (Fig. [I]). 

Phase-diagram of H 0 g: To deduce the structure of 
the h = 0 phase diagram of H ef j (Fig ([ 2 ] A)) in the T- 
•Asing plane (with J xy = 1) for fixed 0(1) values of Jg T 
and he, I start with the known phase diagrams of H xy 
and Hising, and analyze the effects of a non-zero Jg T . 
To this end, recall that Hi s i ng develops long-range or¬ 
der in r for T < T t , with long-distance properties of 
the critical point at T t described by a fixed-point free- 
energy functional F 1 / 2 = J d 2 xHi/ 2 , with central charge 
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FIG. 1: Color-coded symbols on sites give the value of (S~) 
in the presence of ferrimagnetic {9 = 0) or antiferromagnetic 
(9 = three-sublattice order in spin-5 triangular [Kagome] 
lattice easy-axis antiferromagnets. These ordering patterns 
distinguish between the three sublattices of the underlying 
Bravais lattice of sites [up-pointing triangles]. 



FIG. 2: (Color online) A) Predicted structure of the T-Jising 
phase diagram of H e g for h = 0 and fixed J xy and Jg T . Phase 
boundaries of H e g are depicted by colour-coded solid lines, 
while those of Hi S mg and H xy are displayed as dashed lines. B) 
Known T-hz phase diagram of H xy + hz ^]-cos(3 6?) showing 
the three-state Potts line T c {hz). Path L in A) maps to the 
eponymous path in B). 


c = 1/2. Similarly, H xy develops six-fold symmetry¬ 
breaking long-range order in 9 for T < Tg i, which 
melts via an intermediate phase with power-law corre¬ 
lations: (.s „ i/|f| r) ( T ) wither) e (|,|) 
for temperatures T £ (Xgi, Tg 2 ) HE J3] 02]. Long- 
wavelength properties of this power-law ordered phase 
are controlled, in renormalization group (RG) language, 
by a c = 1 line of fixed points [32], with effective free- 
energy Fkt = / d 2 rJ kt, where 


?KT /T 


y^(W ) 2 

4ng y 


(3) 


with g(T) £ (|, |) corresponding to T £ {Tg\,Tg 2 ). This 
fixed-line has power-law correlations (e l ( 6 , b’ ) - 0 (o'))^ ^ 
1 J r v(g) w ith g{g) = g, which render the six-fold 
symmetry-breaking perturbation h§ cos( 6 dp) irrelevant 
for g > 1/9, and vortices in 6 irrelevant for g < 
1/4 |32j. However, the three-fold symmetric perturbation 
/i3cos(30f) is relevant everywhere on this fixed line [52] , 
implying that long-range order sets in at infinitesimal hz 
when T < Tg 2 . In contrast, for fixed T > Tg 2 , long-range 
order sets in via a three-state Potts transition [39] (Fig. [2] 
B) only when a threshold /i 3 c (T) is crossed; this defines a 
three-state Potts critical line T c (hz ) in the (T, / 13 ) phase 
diagram of H xy (Fig. |2] B)). Therefore, our analysis 
splits naturally into two cases, T t < Tg 2 , and T t > Tg 2 , 
and relies crucially on the observation that long-range 
order of 8 in H xy leads to an external magnetic field 
of effective strength h e g = Jg T (cos(39)} acting on r in 
Hising, while long-range order of r in Hi slng perturbs H xy 
by a three-fold symmetric term X/,-?/i 3 e ff cos(36V), with 

^3eff = 

Tt < Te 2 : If Hising is in a short-range correlated 
paramagnetic phase in the entire temperature range 
(Tgi,Tg 2 ), i.e. if T r < Tgi , a non-zero Jg T only renormal¬ 
izes the value of g(T) that controls the power-law corre¬ 


lators of 9 in this regime. And when the temperature 
is lowered below Tgi, long-range order of 9 in H xy gives 
rise to an effective field h e g = Jg T (cos(39)) in Hi smg , con¬ 
verting the Ising transition at T t to a smooth crossover. 
On the other hand, if Tg\ < T Tl long-range order of r 
below T t leads to a three-fold symmetric perturbation 
/i 3 eff of Hxy, which immediately causes H xy to develop 
long-range order in 9 (Fig. [ 2 ] B). 

Thus, when T t < Tg 2 , H 0 g is expected to display a six¬ 
fold symmetry-breaking long-range ordered state for T < 
T c 1 , which undergoes a two-step melting transition via 
an intermediate power-law ordered phase (corresponding 
to T £ (T c i,T c2 )) with an exponent g(T) that increases 
from rj(T c 1 ) = | to g(T c2 ) = The value of T c 1 is set 
(with deviations of order Jg T ) by the larger of Tg\ and 
T t , while that of T c2 is approximately set by Tg 2 . 

Power-law ordered phase: Long-wavelength prop¬ 
erties of H e g in this power-law ordered intermediate 
phase can be described quite generally (for either sign 
of he ) by an effective free-energy density 

TtKt/T = Tkt /T + cg T rp cos(3#p) . (4) 

Although a nonzero cg T leads, upon tracing over r, to 
the six-fold term cos(60^) which is irrelevant all along 
the fixed-line parametrized by cg T = 0 and g{T) £ (|, 4) 
[as in in Eqn. pj], I choose to retain a bare c T g ■=£ 0 ex¬ 
plicitly in Eqn. [4] since this “dangerously irrelevant” cou¬ 
pling controls the long-distance correlations of Tp along 
this fixed-line. Indeed, the nonzero value of cg T in 3 >kt 
causes Tp to inherit the power-law correlations of cos(3 9p) 
for all T £ (T c i, T c2 ): (r r -r 0 ) ~ ( e 3i ( e ~ e G) „ l/ r "( T \ 
Ferromagnetic couplings between the Ising spins are 
not explicitly included in 3 >kt/T since the Ising bond- 
energy Eip p 2 \ = Tp Tf 2 has rapidly decaying correla¬ 
tions ~ 1 /r 369 (r is the distance between 

bonds (fif 2 ) and (f^r/)) that render these couplings ir¬ 
relevant along this fixed line. Just below g = 1/9 {i.e. 
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FIG. 3: L dependence of C t (tl) and C p g{f L ) (p=l,2,3) 
at separation tl = e x on periodic L x L triangular 
lattices, evaluated at the estimated location [//(),/j M ] = 
[1.5570,1.0061] of the multicritical point of H c g with J xy = 
h 6 = 1.0, Jg T = 0.25 (notation as in text). Lines denote fits 
to 1 / L Vt and 1 /L v P e respectively, using 773 s = g T = 0 . 201 , 
77 s = 0.258, and 772 s = 0.353. Cie [C 30 ] is rescaled by a factor 
of 7 [factor of 10] for clarity. 

for T < T c 1 ), the ferromagnetic couplings between the rp, 
and the six-fold anisotropy term cos(6(9p), both become 
relevant. This signals the onset of six-fold symmetry¬ 
breaking long-range order in H e g. 

Singular susceptibility: For r](T) < | in this 
power-law ordered intermediate phase of H e g, the fore¬ 
going implies that power-law correlations of r decay 
slowly enough that they lead to a divergent contribu¬ 
tion Xsing. ~ L 2 ~ 9v to the finite-size susceptibility \l 
of an L x L system at h = 0. This implies Xl{T) = 
Xreg (T) + b(T)L 2 ~ 9 ^ for 7 7 (T) £ (£,§). When an 
external field h is turned on in this regime, it per¬ 
turbs 3>kt with a three-fold symmetric perturbation 
J 0r Xreg/t cos(3dy). This drives H c g to a long-range or¬ 
dered state with correlation length £(h) ~ |/i| -1 /* 3 ( fl ), 
where A 3 (g) = 2 — 9g/2. Beyond this correlation-length 
scale, H b s resembles a three-state Potts model in its or¬ 
dered state [23J. Therefore, for small non-zero h , Xsing. 
will be cut off at length-scales of order this correlation 
length £(h), giving rise to a thermodynamic susceptibility 
that scales as (£( h )) 2 ~ 9ri< - T ) at small h. For the thermo¬ 
dynamic easy-axis susceptibility of the microscopic easy- 
axis antiferromagnet, the foregoing analysis thus predicts 

X(B) ~ \B\ 4-viVn (5) 

at small \B\ for 77 (F) £ (|, |). This prediction identi¬ 
fies an experimentally useful signature of two-step melt¬ 
ing of either type (ferrimagnetic or antiferromagnetic) 
of three-sublattice order in triangular and Kagome lat¬ 
tice easy-axis magnets. In particular, it applies to the 
S = 1 triangular lattice Heisenberg antiferromagnet with 


strong single-ion anisotropy [2311333, and to the triangular 
lattice Ising antiferromagnet with further-neighbour cou¬ 
plings pJ2j. It would therefore be interesting to identify 
quasi two-dimensional magnets in the Ca 3 Co 206 fam¬ 
ily [511 l 5Tj (with an angular momentum J = 1 ion at 
one Co site and a nonmagnetic ion at the other) which 
could provide experimental realizations of the former. It 
would also be interesting to identify new combinations 
of substrate and adsorbate for which monolayer densities 
closer to half-filling (than hitherto achievable [511111 ). cor¬ 
responding to B <C 1 in the latter, could be reached for 
monolayers of adsorbed gases on triangular substrates. 

T r > T 02 : In this case, H c g develops long-range or¬ 
der in r via a transition in the Ising universality class 
at T cI (Fig [2] A), with the value of T cI set by T t (with 
deviations of order Jg T ). For T < T c i, the spontaneous 
magnetization (r) perturbs H xy with the three-fold field 
/i 3 e ff. Lowering the temperature below T c j along path 
L in the phase diagram of H e g (Fig. | A) therefore cor¬ 
responds to moving along the eponymous path L in the 
known [39j| phase diagram of H xy + h 3 JWcos(3 dp) (Fig. [ 2 ] 
B). This key observation immediately leads to two con¬ 
clusions: First, H e g must develop long-range order in 6 
at a lower temperature T c p < T c j via a three-state Potts 
transition (Fig. HA). Second, these Ising and three-state 
Potts transition lines (T c j and T c p) must meet the phase 
boundaries of the power-law ordered phase (T C 2 and T c i) 
at a single multicritical point M (Fig. |A). 

Multicritical point: The fixed point theory Fjg that 
controls long-distance properties of M can be reached 
from the c = 3/2 theory F 1 / 2 + F K t (with g = 1/4) by 
turning on the relevant perturbation Jg T . The c-theorem 
[44] therefore predicts that the central charge of Fm 
obeys cm < | ■ Since Fm must have a relevant direction 
leading from it to the c = 1 theory F T kt, the c-theorem 
also predicts Cm > 1. At M, the correlation func¬ 
tions C T (r) = (r(r)r(0)) and C p g(r) = (e^Me-^ 0 )) 
(p = 1,2,3) are expected to have the long distance forms: 
C T (r) = l/r VT , C p g(r) ~ 1 /r ripS (with 773 g = r] T on sym¬ 
metry grounds). Setting J xy = he = 1.0, Jg T = 0.25 
and parametrizing Ji sing = f xy T ei /T T and T = fif xy T gi , 
with Tgi = 1.04 and T t = 3.6409, I have performed ex¬ 
tensive Monte-Carlo simulations of H 0 g to locate and 
study M. Fig [3] displays power-law fits for the L depen¬ 
dence of C T {fi,) and C p q{tl ) at separation tl = on 
periodic Lx L triangular lattices at my best estimate for 
M, given by ~ [1.5570(8), 1.0061(5)]. Such fits 

yield the following estimates for multicritical exponents: 

773 , = Vt « 0.201(20); 77 , » 0.258(5); 772 g « 0.353(6). (6) 

This set of exponents is clearly different from the well- 
known exponents in the power-law ordered phase [ 77 ^, = 
p 2 77(T)], or on the three-state Potts line [ 772 © = 77 g = 4/15] 
or the Ising line [ 773 , = r] T = 1/4]. 

I close by noting an intriguing possibility: Since three- 
sublattice order melts via a three-state Potts transition 
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in the Kagome Ising antiferromagnet with dipolar inter¬ 
actions [MU], while the analogous short-ranged model 
with nearest and next-nearest neighbour exchange cou¬ 
plings exhibits two-step melting behaviour mmm, 
it appears likely that such multicritical melting could be 
seen in artificial Kagome-ice systems if the strength of the 
first and second-neighbour exchange interactions could 
be increased relative to the long-range dipolar couplings 
(whose values are fixed by magnetostatics). 
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in the Program on Frustrated Magnetism and Quantum 
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was informed by useful discussions with T. Grover and a 
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ner. Participation in this program was made possible by 
partial support from National Science Foundation grant 
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